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ABSTRACT 

The very high energy (VHE; E > 100 GeV) blazar Markarian 501 has a well-studied 
history of extreme spectral variability and is an excellent laboratory for studying the 
physical processes within the jets of active galactic nuclei. However, there are few 
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detailed multiwavelength studies of Markarian 501 during its quiescent state, due to 
its low luminosity. A short-term multiwavelength study of Markarian 501 was coordi- 
nated in March 2009, focusing around a multi-day observation with the Suzaku X-ray 
satellite and including 7-ray data from VERITAS, MAGIC, and the Fermi Gamma-ray 
Space Telescope with the goal of providing a well-sampled multiwavelength baseline 
measurement of Markarian 501 in the quiescent state. The results of these quiescent- 
state observations are compared to the historically extreme outburst of April 16, 1997, 
with the goal of examining variability of the spectral energy distribution between the 
two states. The derived broadband spectral energy distribution shows the character- 
istic double-peaked profile. We find that the X-ray peak shifts by over two orders of 
magnitude in photon energy between the two flux states while the VHE peak varies 
little. The limited shift in the VHE peak can be explained by the transition to the 
Klein-Nishina regime. Synchrotron self-Compton models are matched to the data and 
the implied Klein-Nishina effects are explored. 

Subject headings: BL Lacertae objects: individual (Markarian 501 = VER J1653+397) 
— galaxies: active — gamma rays: observations — radiation mechanisms: non-thermal 
— X-rays: galaxies 



1. Introduction 

Blazars, a subclass of active galactic nuclei (AGN), are the dominant extragalactic source class 
in 7-rays. They have been observed to show rapid variability and non-thermal spectra, presenting a 
broad continuum across nearly the entire electromagnetic spectrum. This implies that the observed 
photons originate within highly relativistic jets oriented very close to the observer's line of sight 



(jUrrv Sz Padovanil 119951 ). This orientation results in Doppler beaming that boosts the intensity 
and frequency of the observed jet emission, often overwhelming all other emission from the source. 
Therefore, blazars make excellent laboratories for studying the physical processes within the jets 
of AGN. They were among the first sources to be detected in the VHE band, and as of this writing 
there are 34 known VHE 7-ray blazarg^. Markarian 501 (Mrk 501; 1H1652 +398), at a redshift of 



z = 0.034, was the second blazar to be detected at VHE (|Quinn et al.lll996l ). 



The spectral energy distribution (SED) of Mrk 501 characteristically shows a double-peaked 
profile. These peaks occur at keV and TeV energies when the SED is plotted in the i>~F v vs v 
representation. This structure is common among all VHE 7-ray blazars, and several models have 
been developed to account for the double-peaked structure. These models uniformly attribute the 
peak at keV energies to synchrotron radiation from relativistic electrons and positrons within the 
blazar jets, but they differ in accounting for the source of the VHE peak. The models are generally 



http: / /tevcat . uchicago.edu / 
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divided into two classes: leptonic and hadronic, named for their attributed source for the VHE 
peak. The leptonic models advocate inverse-Compton scatt ering to VHE of either the synchrotron 



photons from within the jet or an external ph oton field (e.g.. lMarscher Gearlll985l ; lMaraschi et al 



1992 ; iDermer et al.lll992l : ISikora et al,lll994l ). The hadronic models, however, account for the VHE 
emission by tt° decay or by ir^ decay with subsequent synchrotron and/or Compt on emission from 



decay products, or by synchrotron radiati on from ultra-relativistic protons (e.g.. iMannheim 



1993 



Aharonian 



2000 



Pohl fc Schlickeiserll2000h . 



Observationally, Mrk 501 has been known to undergo both major outbursts on long time 
scales and rapid flar e s on short time s cales, most prom i nently in the keV a nd VHE range (e.g., 
Catanese et all Il997l : IPian et al.1 Il998l : IXue & Cuil hood : lAlbert et all l2007ah . During these out- 
bursts, both of the SED peaks have been observed to shift towards higher energies. During the 
most extreme cases, the keV peak has been observed above 200 keV, well above typical values be- 
low 1 keV. Historically, the SED has been measured in the VHE band primarily during outbursts, 
due to the lower sensitivity of previous generations of instruments. A previous study examined 
the quiet state, but it was performed before the launch of the Fermi Large Area Telescope (Fermi - 
LAT), which provides coverage in the range between keV and VHE energies (jAnderhub et al.l l2009). 
This work attempts to provide state-of-the-art short-term multiwavelength measurements of the 
quiescent state of Mrk 501, with broad spectral coverage in the critical keV and VHE bands as 
well as coverage with the Fermi-L AT . These measurements are then compared to observations by 
BeppoSAX in X-rays and by the Whipple 10m and CAT Cerenkov telescopes in VHE 7-rays during 
the 1997 extreme outburst. This outburst has been well studied using multiple inst ruments (e.g., 
Aharonian et al.ll200ll : ICatanese et al J 1 199 71 : 1 Aharonian et al.ll200ll : IPian et al.lll998l ) and provides 
a good comparison to the quiet state observed in 2009. 



2. Data Selection and Analysis 
2.1. Suzaku: X-ray 



The Suzaku X-ray observatory, a collaborative project between institutions in the United States 
and Japan, is an excellent tool for studying the broadband SED of sources such as Mrk 501 due to 
its broad energy range (0.2 - 600 keV). Suzaku has two operating in struments for studying X-ray 
emissio n: the X-ray imaging spectrometer (XIS: lKovama et al.ll2007l ) and the hard X-ray detector 
(HXD; iTakahashi et"al"1l2007h . 



The XIS instrument consists of four X-ray telescopes with imaging CCD cameras at their focal 
planes. Three of the CCDs are front-illuminated, and one CCD is back-illuminated to provide 
extended sensitivity to soft X-rays. The combined energy range of these CCDs is 0.2 - 12.0 keV. 
The HXD is a non-imaging instrument that expands the energy sensitivity of Suzaku to the 10 - 
600 keV band by using two types of detectors. Silicon PIN diodes provide sensitivity in the range 
10 - 70 keV, and gadolinium silicate [GSO; Gd2SiOs(Ce)] scintillators placed behind the PIN diodes 
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provide sensitivity in the range 40 - 600 keV. One limitation of Suzaku with respect to attaining 
simultaneous multiwavelength observations is that it resides in a low earth orbit, so observations 
are subject to frequent earth occultations, limiting continuous temporal coverage. 

Observations of Mrk 501 were carried out with the Suzaku X-ray satellite from 2009-03-23 UT 
18:39 to 2009-03-25 UT 07:59 (sequence number 703046010). After good time interval selection 
and dead-time correction, a total of approximately 72 ks of live time remained for the XIS and 61 
ks for the HXD. To reduce the occurrence of photon pile- up in the XIS, a common problem when 
dealing with bright X-ray point sources, the XIS observations were carried out in 1/4 window mode 
and Suzaku was operated with a pointing nominally centered on the HXD field of view (FoV). In 
addition, HXD nominal pointing increases the effective area of the HXD instrument. The resulting 
XIS count rates are well below the threshold where pile-up becomes a significant issue. Even with 
the HXD nominal pointing, there was no significant detection in the GSO scintillators, so the 
maximum useful energy for these observations using only the HXD/PIN detector was around 70 



The raw XIS data were reprocessed following the guidelines from the Suzaku teano using 
Suzaku FTOOLS version 12.0. The events were graded and only those with ASCA grades 0, 2, 
3, 4 and 6 were included in further analysis. In addition, the data were screened to exclude times 
when the Suzaku spacecraft was passing through the South Atlantic Anomaly (SAA), including 
a 436 s post-SAA buffer to allow the instrument backgrounds to settle. Additional cuts were 
made requiring the data used to be taken at least 5° in elevation above the Earth's limb and at 
least 20° above the Earth's daytime limb. These cuts help reduce contamination from the Earth's 
atmosphere. 

Source events were extracted from a circular region with a radius of 1.92' centered on Mrk 
501. Background events were extracted from a circular region with the same radius located far 
from the source in each XIS detector. The redistribution matrix file (RMF) and ancillary response 
files (ARF) were generated for each XIS detector using the analysis tools included in the Suzaku 
analysis package. 

The HXD data were also reprocessed from the raw files, following the guidelines provided by 
the Suzaku teamjf] The data were filtered with similar criteria to the XIS data with a slightly longer 
settling time of 500 s required following passage through the SAA. In addition, the only elevation 
requirement was for data to be taken at least 5° above the Earth's limb. 

Because the HXD is a non-imaging instrument, the cosmic X-ray background (CXB) and 
instrumental non- X-ray background (NXB) must be accounted for in different ways than for the 
XIS. The time- variable non-X-ray background has been modeled and provided by the HXD team. 
The cosmic X-ray background is not included in this model and must be accounted for separately 



2 http://heas arc. gsfc.nasa.gov/docs/suzaku/analysis/abc/node9. html 
3 http://heas arc. gsfc.nasa.gov/docs/suzaku/analysis/hxd_repro. html 



keV. 
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as the CXB flux is about 5% of the background for PIN observations. The model suggested by 
the HXD team is based on results from HEAO-1 and modified for inpu t to XSPEC . The imported 
model reproduces the PIN background with an HXD nominal pointing (|Boldtlll987l ) . The model is 
given by: 



dF { E \ / jtj 

— = 8.0 X 10- 4 x exp (4=^) Photons cm" 2 s 1 KOV" 1 koV" 1 . (I) 

The CXB background model was added to the NXB background model and used to sub- 
tract background events in XSPEC during spectral analysis of the PIN data. The response file 
ae_hxd_pinhxnome5_20080716 . rsp provided by the HXD team was used for the PIN spectral anal- 
ysis. 



2.2. VERITAS: VHE 7 -ray 



The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is an array of 
four imaging atmospheric-Cerenkov telescopes (IACTs) located at the Fred Lawrence Whipple 
Observatory in sou thern Arizona (31°40'N, 110°57'W) at an altitude of 1.28 km above sea level 
([Holder et al.ll2008l ). These telescopes capture the Cerenkov light emitted by extensive air show- 
ers that are initiated in the upper atmosphere by 7-rays and cosmic rays. The telescopes have 
12m-diameter tessellated reflectors directing light onto imaging cameras, each consisting of 499 
photomultiplier tubes (PMTs), located in the focal plane of each telescope. The cameras have a 
FoV of 3.5° with an angular size per pixel of 0.15°. VERITAS is sensitive in the range of approxi- 
mately 100 GeV - 30 TeV with an energy resolution of between 15 and 20%. Images of the extensive 
air showers captured by each telescope can be combined to reconstruct the direction and energy of 
the shower initiator. The telescopes typically operate in "wobble" mode, where the location of the 
object to be observed is offset fr om the center of th e FoV by 0.5°, allowing for simultaneous source 
and background measurements (jFomin et al.lll994i ). The offset direction alternates between north, 
south, east and west for each data run (typically lasting 20 minutes) to reduce systematic effects 
from the offset. 



The VERITAS data were analyzed using the standard analysis package (jCogan et al.l 120071 ). 
Images of the Cerenkov light produced by the extensive air showers are parameterized, and quality 
cuts are placed on the data to ensure that the shower reconstruction procedure can be performed 
reliably. Separation of showers initiated by 7-ray photons from hadronic showers initiated by cosmic 
rays is performed with the aid of lookup tables. The lookup tables are calculated for each telescope 
with several values of zenith angle, azimuth angle, and the night sky background noise level. The 
observed parameters length and width are scaled by the corresponding parameters in the lookup 
tables and averag ed among the teles copes to find the "mean scaled parameters" in the method 
first suggested by iDaum et al.1 (| 19971 ). Selection criteria for 7/hadron separation (optimized with 
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data taken on the Crab Nebula) are placed on the mean scaled parameters of the shower images to 
exclude most hadronic showers while still retaining a large proportion of 7-ray initiated showers. 

To perform a background subtraction of the surviving cosmic-ray event s, an estimation of these 
backg round counts is made using the "reflected-region" background model (IBerge. Funk, Hinton 



20071 ). Events within a squared angular distance of 6 2 < 0.0169 of the anticipated source location 



1983 ). More 



Acciari et al 



are considered ON events. Background measurements (OFF events) are taken from regions of the 
same size and at the same angular distance from the center of the FoV. For this analysis, a minimum 
of eight background regions were used. The excess number of events from the anticipated source 
location is found by subtracting the number of OFF events (scaled by the relative exposu re) from 
the ON events. Statistical significances are calculated using Eqn. 17 of iLi &: Ma 
detail s about VERITAS, the calibration and analysis techniques can be found in 

liooah . 

On March 24, 2009 (MJD 54914), VERITAS was used to take data on Mrk 501 while operating 
in wobble mode, with an offset angle of 0.5°. VERITAS was used to observe Mrk 501 for 2.6 hours 
from 9:11 UT to 11:58 UT using a series of 20-minute runs, with a live-time of 2.4 hours and zenith 
angles ranging from 35° to 9°. Additional observations totaling 2.9 hours were made on March 23 
and 25 but were discarded due to poor weather conditions. The data were processed using standard 
VERITAS analysis cuts, finding 431 ON events and 1810 OFF events with a maximum ratio between 
ON and OFF exposure of 0.125. The total significance of the excess for this observation was 12.7 
standard deviations (a) with an average significance per run of 3.8 a. 



2.3. MAGIC: VHE 7-ray 

MAGIC is a system of two IACTs for VHE 7-ray astronomy located on the Canary island of 
La Palma (2.2 km above sea level, 28°45'N, 17°54'W). At the time of the observation, MAGIC-II, 
the second telescope of the array, was still in its commissioning phase so Mrk501 was observed by 
MAGIC-I in stand-alone mode. 

MAGIC-I has a 17m-diameter tessellated reflector dish and a hexagonal camera with a FoV 
of 3.5° mean angular diameter. The camera comprises 576 high-sensitivity PMTs of different 
diameters (0.1° in the inner portion of the camera and 0.2° in the outer portion). The trigger 
threshold of MAGIC at low zenith angles is between 50 and 60 GeV with the standard trigger 
configuration used in this campaign. The accessible energy range extends up to tens of TeV with 
a typical energy resolution of 20% to 30%, depending on zenith ang le and energy . Furth er details, 



telescope parameters, and performance information can be found in lAlbert et al.1 (j2008al ) 



The MAGIC data for this campaign were recorded between 3:50 UT and 5:10 UT on the night 
of March 23, 2009 (MJD 54913). The observations were performed in wobble mode. The zenith 
angle of the observations ranged from 15° to 30°. The recorded data fulfil the standard quality 
requirements, so none of the runs were removed from the data sample. The total live time for this 
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observation is 1.3 h. 



The data were analyzed using the standard MAGIC calibration and analy sis ([Albert et al 



2008al ). The analysi s is based on image parameters (jHillas 



1985 



Aliu et al.ll2009l ) and the random 



forest (RF) method (lAlbert et al.ll2008bl ). which are used to define the so-called hadronness of each 
event. The cut in hadronness for background rejection was chosen using Monte Carlo data, setting 
a cut efficiency for 7-rays of 70%. Three background (OFF) sky regions are used for the residual 
background estimation, giving a ratio between ON and OFF exposure of 0.333. The OFF regions 
are chosen to be at the same angular distance from the camera center as the ON region. The excess 
of 7-ray events from the source was extracted from the distribution of the ALPHA parameter, 
which is related to the image orientation on the camera (jHillas! Il985l ). The final cut on ALPHA 
was also chosen using Monte Carlo 7-ray events, again imposing a cut efficiency of 70%. With these 
less restrictive cuts, suitable for the spectrum calculation, the significance of the excess is 2.2 a 
with 9527 ON events and 27846 OFF events, whereas applying stricter cuts in the optimal energy 
range for signal detection leads to a significance of 9.7 a with 123 ON events and 97 OFF events. 

The energies of the VHE 7-ray candidates were also estimated using the RF method. A cut 
was applied on the SIZE parameter, requiring a minimum of 100 photoelectrons in the shower 
images. Because the SIZE parameter is related to the energy of the initiator of the shower, this 
cut results in an energy threshold for the analysis (defined as the peak position of the distribution 
of the reconstructed energy parameter for simulated 7-ray events after the cut) of ~ 100 GeV. 



2.4. Fermi-LAT: HE 7-ray 

The Fermi Large Area Telescope (LAT) is a satellite instrument to perform 7-ray astronomy 
from 20MeV to several hundred GeV. The instrument is an array of 4 x 4 identical towers, each 
one consisting of a tracker (where the photons are pair-converted) and a calorimeter (where the 
energies of the pair-converted photons are measured). The entire instrument is covered with an 
anticoincidence detector to reject charged-particle background. The LAT has a peak effective area 
of 0.8 m 2 for 1 GeV photons, energy resolution typically better than 10%, and a FoV of about 2.4 sr, 
with angular resolution (68% containm ent angle) better tha n 1° for energies above 1 GeV. Further 



details about the LAT can be found in lAtwood et al.l (120091 ) 



The analysis was performed with the ScienceTools software package version v9rl5p2, which is 
available from Fermi Science Support Center (FSSC). Only events having the highest probabilities 
of being photons, those in the "diffuse" class, were used. The Pass 6v3 diffuse-class instrument 
response functions were used for the analysis. The LAT data were extracted from a circular region 
of 10° radius centered at the location of Mrk501. The spectral fit was performed using photon 
energies greater than 0.3 GeV in order to minimize systematic effects, and a cut on the zenith 
angle (> 105°) was also applied to reduce contamination from Earth-albedo 7-rays, which are 
produced by cosmic rays interacting with the upper atmosphere. 
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The background model used to extract the 7-ray signal includes a Galactic diffuse emission 
component and an isotropic component. The model that we adopted for the Galactic component is 
glLiem _v02.fitB The isotropic component, which is the sum of the extragalactic diffuse emission and 
the residual charged-particle background, is parametrized here with a single power-law function. 
Owing to the relatively small size of the region analyzed (radius 10°) and the hardness of the 
spectrum of Mrk501, the high-energy structure in the standard tabulated isotropic background 
spectrum, isotropic_iem_v02.txt, does not dominate the total counts at high energies. In addition, 
we find that for this region a power-law approximation to the isotropic background results in 
somewhat smaller residuals for the model overall, possibly because the isotropic term, with a free 
spectral index, compensates for an inaccuracy in the model for the Galactic diffuse emission, which 
is also approximately isotropic at the high Galactic latitude of Mrk501 (b ~ 39°). In any case, the 
resulting spectral fits for Mrk501 are not significantly different if isotropic_iem_v02.txt is used for the 
analysis. To further reduce systematic uncertainties in the analysis, the photon index of the isotropic 
component and the normalization of both components in the background model were allowed to 
vary freely during the spectral point fitting. In a ddition, the model also includes five nearby 
sources from the 1FGL catalog dAbdo et alJ boiOaJ ): 1FGL J1724.0+4002, 1FGL J1642.5+3947, 
1FGL J1635.0+3808, 1FGL J1734.4+3859, 1FGL J1709.6+4320. The spectra of those sources 
were also parameterized by power-law functions, whose photon index values were fixed to the 
values from the 1FGL catalogue, and only the normalization factors for the single sources were left 
as free parameters. The spectral analysis was performed with the post-la unch instrument-res ponse 
functions P6_V3_DIFFUSE using an unbinned maximum likelihood method ([Mattox et all 19961 ). The 
Fermi collaboration estimates the current systematic uncertainties on the flux as 10% at 0.1 GeV, 
5% at 560 MeV and 20% at 10 GeV and above@. 



3. Results 
3.1. Light Curves 

The time coverage of the data sets used for this campaign is illustrated in Figure [TJ Light curves 
were produced in the passband of each instrument involved in the campaign to understand the flux 
levels and variability of Mrk501 during these observations; however, the Fermi-LAT detection was 
not significant enough to include a light curve. The resulting multiwavelength light curves are 
shown in Figure [2j 

Light curves produced from the XIS and HXD/PIN Suzaku data, with 5 ks time bins, are 
shown in the top two panels of Figure [2j Individual light curves for each of the XIS detectors were 
created and then summed. The best fit of a constant for the summed XIS (0.6 - 10.0 keV) light 



|http://fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 



5 See http : //f ermi . gsf c . nasa. gov/ ssc/data/analysis/LAT_caveats .html 
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curve is 16.14 counts s _1 with a \ 2 value of 4207 for 26 degrees of freedom. The poor fit of a 
constant to the data corresponds to the moderate variability that is visible by inspection in the 
XIS count rate. A constant value was also fitted to the HXD/PIN light curve, yielding a value of 
0.11 counts s _1 with a x 2 value of 43.7 for 26 degrees of freedom, again indicating variability. 

In the VHE 7-ray band, the source appears to be in a "quiescent" st ate when compared 
to contemporaneous long-term stud ies of the source's broadband behavior (jAbdo et all l2010bl ; 
Huang et allbood : IPichel et al.lbood 'l. The light curves produced from the VERITAS and MAGIC 
data used 20 minute time bins and are plotted in the bottom panel of Figure [2j To produce light 
curves, the spectral slope from the corresponding data is assumed (see Section f3.2|) . and the effec- 
tive collection areas were used with the excess counts to determine the spectral normalization. An 
energy threshold of 300 GeV was then used to calculate the integral flux for each bin. A constant 
flux value was fitted to the VERITAS and MAGIC data separately, yielding a best fit value of 
(2.8 ± 0.5) x lO" 11 photons cm' 2 s" 1 for VERITAS and (3.9 ± 0.8) x lO" 11 photons cm" 2 s" 1 for 
MAGIC. The VERITAS data show no evidence for any variability, with a x 2 value of 3.6 for 7 
degrees of freedom for the constant fit. Likewise, the constant fit to the MAGIC data yields a x 2 
value of 3.2 for 3 degrees of freedom. The individual constant fits to the VERITAS and MAGIC 
data differ slightly but show that the flux measurements from the two instruments are reasonably 
compatible, with no large-scale variability. 



3.2. Spectral Analysis and Modeling 



The spectral results for the observations taken during this campaign when Mrk 501 was in 
the quiescent or low-flux state are shown in Figu res El IH an d JH F o r comparison, t h ey ar e also 
shown in Figure El along w ith archival data from iPian et all (| 19981 ) , ICatanese et al.1 (119971 ) , and 
Djannati-Atai et al.1 (|1999l ) taken when Mrk 501 was in a high-flux state. 



The data from all Suzaku instruments were fitted jointly in XSPEC 12 (|Arnaudl ll996), A 
cross-normalization parameter was introduced to account for the calibration differences between the 
individual XIS detectors and the HXD/PIN. A broken power-law function modified by interstellar 



absorption with a fixed column density N# = 1.73 x 10 cm (jStark et al.lll997l ) was fitted to 



the data, yielding power-law indices Ti = 2.133 ± 0.003 s t a t and T2 = 2.375 ± 0.009 s t a t, and a 
break energy of Eb rea k = 3.21 ± 0.07 s t a t keV. After the fitting procedure, the cross-normalization 
factors for the XIS1 and XIS3 detectors relative to the XISO detector are c ompatible with previous 
calibration measurements of the Crab Nebula flux (jSerlemitsos et al.ll2007l ). A fixed normalization 
factor of 1.18 was introduce d between the HXD /PIN and XISO as specified by the Suzaku team 
for HXD nominal pointings ( Maeda et al.ll2008l ). The resulting fit, with x 2 /d.o.f.= 1.33 for 1511 
degrees of freedom, is plotted in Figure El After an acceptable fit was found, the result was unfolded 
through the instrument response and de-absorbed to derive the intrinsic X-ray spectrum shown in 
Figure El 
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The Fermi-h AT spectral measurement shown in Figures [5] and [6] is derived using a time interval 
of 7 days centered on the Suzaku observing time window - from MJD 54911 up to MJD 54918. 
Mrk501 is detected with a Test Statistic (TS) value of 46 (~ 7a). The spectrum can be fitted with 
a single power-law function with photon flux F(> 0.3 GeV) = (1.7 ± 0.9 s t a t) x 10 -8 ph cm _2 s _1 
and differential photon spectral index Flat = 1.7 ± 0.3 s t a t. The number of photons (predicted by 
the model) for the entire spectrum is 12, out of which 9 are in the energy range 0.3-3 GeV, 2 in the 
energy range 3-30 GeV and 1 in the energy range 30-300 GeV. In order to increase the simultaneity 
of the reported SED, we also performed the analysis of the LAT data over a period of 3 days 
centered at the Suzaku observations, namely the time interval from MJD 54913 to MJD 54916. 
For such a short time interval the detection has a TS value of 10 (~ 3a), and significantly greater 
uncertainties in the power-law parameters: F(> 0.3 GeV) = (1.8 ± l-3 s tat) x 10 _8 ph cm -2 s~ , 
F LAT = 1.9 ±0.5 

stat- Since the two spectral shapes are compatible within errors and there is no 
indication of flux variability during this time, the spectrum derived with the 7-day time interval is 
used in the SED plot in Figures H] and [6l 

For the VERITAS spectral analysis, the spectrum was unfolded using effective collection areas 
for the array that were calculated using Monte Carlo simulations of extensive air showers passed 
through the analysis chain with a configuration corresponding to the data-taking conditions and 
the 7/hadron selection cuts. The unfolding process takes into account the energy resolution of the 
array and potential energy bias. The effective areas are then used for each reconstructed event to 
calculate the differential photon flux. The VHE 7-ray differential spectral points were fit with a 
simple power law of the form: 



<1E " ' 1 X \lTeV I Photon cm~ z s~ L TeV" 1 . (2) 

For the VERITAS analysis, best-fit parameters of Fy^E = 2.72 ± 0.15 sta j ± 0.l sys and F = 
5.78 ± 0.83 s tat ± l.W sys were found. 

For the MAGIC analysis, the d erived spectrum wa s unfolded to correct for the effects of 



the energy resolution of the detector (j Albert et alj|2007bl ) and of possible bias. Four statistically 



significant spectral bins were obtained. The spectrum is compatible with a simple power-law 
function with T VHE = 2.67±0.21 stat ±0.20 sys and F = 8.34±1.53 stat ±2.50 sys within the MAGIC 
energy range (see Figure [5]) . 

A si mple synchrotron self-Com pton (SSC) modej^is used to match the results with an emission 



scenario (jKrawczvnski et al.l 12004 ) . This model assumes a spherical emission region of radius R 



with the emission Doppler-boosted by the factor: 

D = [T Lorentz (l- PcosO)]- 1 , (3) 



The SSC code is available at http://jelley.wustl.edu/multiwave/spec trum/7codel 
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where the emission region is moving with Lorentz factor T^ oren t z towards the observer, f} gives 
the bulk velocity of the plasma in terms of the speed of light, and 9 is the angle between the 
jet axis and the line of sight in the observer frame. The emission region contains an isotropic 
non-thermal electron population and randomly oriented magnetic field, B. The electron spectrum 
is assumed to be in the form of a broken power-law function with index p\ between 7 m j n and 7^ 
and index p2 bet ween 7h and 7mn.x- Th e m odel accounts for 7-ray absorption by the extragalactic 



background light (jFranceschini et al.ll2008l ). After computing the synchrotron and inverse-Compton 
emission, the model corrects the photon spectrum that is produced for synchrotron-self-absorption 
and internal absorption due to pair-production processes. The model does not evolve the electron 
spectrum self-consistently. 

In general, it is difficult to fully constrain the parameters in SSC models. This requires both the 
flux and energies of the two (synchrotron and SSC) peaks t o be fully determined b y observations 



as well as an independent estimate of the Doppler factor (|Tavecchio et al.lll998l ). In addition, 
measurements of the synchrotron spectrum are used to help determine the indices of the modeled 
electron spectrum (pi and P2). The SEDs from both sets of observations presented here constrain 
the synchrotron and SSC peak locations well (particularly for the low-flux state). In addition, p\ 
for the 1997 high-flux state and P2 for the 2009 low flux-state are well constrained by the data, 
while the other electron spectral indices are not as restricted. The parameters for the models used 
in this work for the 1997 high-flux state and the 2009 low-flux state are given in Table [T] along with 
model parameters from other works for comparison. Some of the model parameters were matched 
to previous models as discussed below. 



4. Discussion 



This data set provides a high-quality sampling of the broadband SED of Mrk 501 in the 
quiescent state, and it allows comparisons to be made over a broad energy range with the extreme 
outburst observed in 1997. In choosing the values for the SED model parameters, we attempted 
to be consistent with previous work in choosing 7 m j n , j ma x, B, R, and D while still matching the 
data (see Table [1]). In addition, when applying the models to the data from the two states, we 
attempted to lim it the differences to the elect ron spectral indices and break energy (75), similar to 
Pian et al.1 (|1998l ) and lAnderhub et al.l (120091 ) . with a small shift in B included to provi de a better 
match to the data. The SSC model for the 1997 outburst was matched to the data fromlPian et al.l 



(jl998l ) and then the shape compared to the data from iDjannati-Ata'i et al.1 (|1999l ). The resulting 
SSC models can reproduce the measured spectra from the keV to the VHE range for both the 
quiescent state and the 1997 outburst. 

The successful match of the SED models to the data from both the outburst and the quiescent 
state primarily by a modification of the injected electron spectrum implies that this may be the 
primary explanation for the dramatic shift in the peak frequency. Specifically, 7^ and the power-law 
spectral index above the break appear to drive the changes. In addition, there is a difference in 
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electron densities, with a density of 2.7 x 10 4 cm~ 2 for the 1997 data and 2.2 x 10 3 cm -2 for the 
2009 data. The SSC model matches indicate that during the high-flux state the break energy is 
shifted an order of magnitude higher and that past the break energy the spectrum is significantly 
harder than in the low-flux state. It should be noted that the match of these models to the data 
is not unique due to the number of free parameters. Although we did not fully explore the multi- 
dimensional parameter space of the model, which is not the main objective of this work, the general 
conclusion on the shift of the SED peaks seems to be robust. 

The X-ray data do not allow us to measure the location of the low-energy SED peaks directly, 
but they can be used to place limits on the peak energies. From inspection of the X-ray data points 
alone, the keV peak should be near 230 keV (5.5 x 10 19 Hz) during the 1997 outburst and near 0.6 
keV (1.5 x 10 17 Hz) during the quiescent state. This means that with a change in flux of around 
one order of magnitude between the quiescent state and the 1997 outburst, the keV peak of the 
SED shifts in frequency by more than two orders of magnitude. 

In contrast to this large shift, the VHE peak of the SED does not seem to shift as dramatically 
in location. The SED peak at VHE is better constrained with the 2009 data than with the 1997 
data, but the SED models applied to the data imply that even with a dramatic shift in flux at 
VHE, the peak energy in the VHE range is stable compared to that in the keV range. This may be 
due to Klein-Nishina (KN) effects which r educe the c r oss-se ction for scattering when dealing with 



electrons of high energy, as also noted by iPian et al.l (| 19981 ). To examine the possible KN effects 



which become important above hv ~ m e c 2 in the electron rest frame, we examined the energies of 
typical electrons from the model distributions. 

From the 1997 SED model, the peak attributed to synchrotron radiation is located at 3.6 x 10 19 
Hz. The peak energy of the emi ssion in the jet frame , takin g the Doppler factor of D = 20 into 
account, is E pea k ~ 7.4 keV. From Rybicki fc Lightman (jl979l ). it can be shown that for synchrotron 
radiation: 

7e ~ 3irqBh ' 1 ' 

where q is charge, j e is the characteristic Lorentz factor of the electrons contributing to the bulk of 
the X-ray emission and E pea k is the peak energy of the synchrotron model. This yields a value of 
7 e = 1.5xl0 6 (£/0.23 G)- 1/2 (D/20r 1/2 for the 1997 outburst. For the 2009 model, the synchrotron 
emission model peaks at 2.24 x 10 17 Hz, giving a peak energy in the jet frame of E pea k ~ 46 eV. 
This yields a value of <y e = 1.0 x 10 5 (£/0.34 G)" 1 / 2 (I?/20)" 1 / 2 for the 2009 measurements. 

Using these electron energies, we can examine the importance of KN effects during the course 
of the two measurements. For the 1997 outburst, 7 e ^ P eafc ~ 11 x 10 4 MeV, well above the electron 
rest mass energy of 0.511 MeV, placing this model scenario in the extreme KN range. For the 2009 
quiescent state model, ^ehfpeak ~ 4.6 MeV, moderately above the KN limit. The extreme energies 
involved during the 1997 outburst may indicate that the emission is piling up in the VHE range 
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due to the reduced KN cross-section. 

The SED models also allow us to examine the characteristic cooling timescales of the syn- 
chrotron and inverse-Compton processes that are assumed to be t he sources of the keV and V HE 



emission, respectively. The synchrotron cooling time is shown by iRybicki &: Lightmanl (| 19791 ) to 
be: 

^ 67rm e c 

where o~t is the Thompson cross-section and 7 e is the characteristic Lorentz factor of the electrons 
contributing to the bulk of the X-ray emission. The inverse-Compton cooling times cale was calcu 



lated for the models, taking into account KN effects using the method discussed in lBottcher et al 



(|1997I ). As discussed above, models for both sets of observations indicate that KN effects are 
important, so these effects were carefully included in the calculation. Using the parameters for 
the 1997 model, the synchrotron cooling timescale is comparable to the SSC cooling timescale 
{T syn ~ 9.8 x 10 3 s, T S sc « 3.7 x 10 3 s) for electrons with 7e = 1.5 x 10 6 (£/0.23 G) - 1 / 2 (13/20) - 1 / 2 . 
However, for the quiescent state model from 2009, the timescale for inverse-Compton cooling is 
two orders of magnitude shorter than the synchrotron cooling timescale (r syn ~ 6.7 x 10 4 s, 
tssc ~ 4.6 x 10 2 s) for electrons with j e = 1.0 x 10 5 (5/0.34 G)" 1 / 2 (L>/20)" 1 / 2 , indicating that in 
this model scenario, radiative cooling is dominated by the SSC process for the values of 7 e calcu- 
lated above. This seems to indicate that the SSC peak energy is determined by the transition into 
the KN regime rather than the maximum or peak electron energies. Beyond this transition energy, 
the ratio of inverse-Compton to synchrotron cooling timescales should decrease dramatically. An 
estimate for the KN transition can be found using the peak of the synchrotron spectrum. For 
example, from the 2009 data, this peak resides at 46 eV in the jet frame. This corresponds to 
a value of e' = E'/(m e c 2 ) = 9 x 10 -5 . The transition to the KN regime is then at the observed 
energy of Ekn ~ m e c 2 D(l/e') ~ 110 GeV, very near the observed SSC peak. In addition, all of 
the timescales calculated were comparable to or less than the light-crossing time for the modeled 
emission regions, indicating that light travel time effects may dominate the light curve profiles. 

Further intense multiwavelength observations of the source, in conjunction with long-term 
monitoring campaigns, will continue to be important to understand the broadband behavior of 
Markarian 501 and to shed light on the blazar phenomenon in general. 
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Observations 
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Fig. 1. — Overall time coverage of Mrk 501 by the instruments involved in the multiwavelength 
campaign. Gaps due to orbital constraints are not indicated in the Fermi and Suzaku plots. A 
color version of this figure is available online. 
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Fig. 2. — Light curves of Mrk 501 in the energy ranges used for Suzaku/XIS (0.6 keV < E < 10.0 
keV), Suzaku/HXD (10 keV < E < 70 keV), MAGIC and VERITAS (E > 300 GeV). Error bars 
show lcr statistical errors. The constant fits to the Suzaku data are shown. See discussion in Sec. 
13.11 A color version of this figure is available online. 



Table 1: Comparisons of SED models for Mrk501 in different flux states from previous work. 
The model s from this work are shown in the last two rows. The 1997 high-flux-state data from 



Pian et al.l (|1998l ) were re-matched for this work with the same SSC code that we used to match 
the data from the 2009 low- flux state (See Sec. 13.21 for a description of the model parameters). 
The models are plotted in Figure [6l 
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Fig. 3. — Observed X-ray spectra of Mrk 501 made from Suzaku XIS and HXD/PIN data separately. 
A joint fit to a broken power- law model was performed for the three XIS CCDs and the HXD/PIN. 
The solid histogram shows the best fits to the data (see Sec. 13. 2\\ . The error bars represent la 
statistical errors. The bottom panel shows the residuals of the fit. 




Fig. 4. — Plot of the Fermi spectral results for the 7 days centered on the Suzaku observing time 
window (MJD 54911 - MJD 54918). The solid line depicts the result of the unbinned likelihood 
power-law fit, and the contour shows the 68% uncertainty of the fit. The data points show the 
energy fluxes computed in differential energy ranges, with a 95% confidence level upper limit for 
the highest energy range. The Fermi analysis was performed in the energy range 0.3-400 GeV. The 
vertical bars show la statistical errors, and the horizontal bars show the energy range for each 
point. A color version of this figure is available online. 
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Fig. 5. — VHE spectrum from VERITAS (diamonds) and MAGIC (triangles), including the best-fit 
power-law models. The VERITAS fit is shown with a solid line, and the MAGIC fit is shown with a 
dashed line. The slopes of the power-law models agree well (T Veritas = 2.72±0.15 S f a t ±0.10 SJ/S and 
T magic = 2.67±0.21<rf a ( ±0.20 sys ). The slight normalization offset corresponds to the moderately 
higher VHE flux levels during the MAGIC observation. The vertical bars show la statistical errors. 
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Fig. 6. — Spectral energy distributions of Mrk 501 in the quiescent state (in filled points) and 
during the 1997 extreme outburst (in open points). The former is made from data obtained in the 
present campaign with Suzaku (filled circles), Fermi (fi lled squares ) , VE RIT AS (filled diamonds). 



and M AGIC (filled triangles). The latter was taken from lPian et al.l (| 19981 ) and lDiannati-Atai' et al 



(jl999l ). with X-ray data from BeppoSAX (open circles), GeV 7-ray upper limit from EGRET (open 
square), and VHE 7-ray data from the Whipple 10m (open diamond) and the CAT Cerenkov 
telescope (open triangles). Note that the Whipple 10m spectral point overlaps with the lowest- 
energy CAT spectral point. The vertical bars for all data show la statistical errors and the upper 
limits are at the 95% confidence level. Corresponding SED models matched to the data for each year 
are shown. The solid lines show the SSC match; dashed lines take EBL absorption into account. 
The VHE data points are not EBL-corrected. A color version of this figure is available online. 



